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ABSTRACT 

y 

Solar flares produce increases in coherent-scatter power from the 
mesosphere due to the increase in free electrons produced by X-ray photo- 
ionization. Thirteen such power enhancements have been observed at Urbana. 
When such an enhancement occurs at an altitude containing a turbulent layer 
with constant strength, we may estimate the relative enhancement of electron 
density from the enhancement in power. Such estimates of enhanced electron 
density are compared with estimates of the X— ray photoionization at that 
altitude, deduced from geostationary satellite measurements* It is found that 
possible types ion-chemical reaction scheme may be distinguished, and the non- 
flare ion-pair production function may be estimated. The type of ion-chemical 
scheme and the nonflare ion-production function are shown to depend on the 
solar zenith angle. 
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INTRODUCTION 

It has, of course, been known for some time that solar flares produce an 
increase in ionization in the D region. This flare-time enhancement in 
ionization results in an increase in coherent scattered power which has been 
assumed to account for an exceptionally full set of good mesospheric velocity 
measurements during the event of April 11, 1978, 0800 CST at Urbana (MILLER et 
al., 1978), and also measured directly for the event of January 5, 1981, at 
1218 AST at Arecibo (ROTTGER. 1983). 

However, the processes linking the solar X-ray enhancement to the co- 
herent-scatter power increase involve many unknowns. Photoionization by X-rays 
may be considered as the driving function of a set of ion-chemical reactions 
which finally determines the electron-density profile; this profile must then 
be advected by turbulence to produce the scattered power. Many details of 
these processes can only be deduced indirectly. 

This paper describes how models of these processes may be constructed 
which account for seme of the features of the power enhancements observed 
during solar flares. Early work along these lines may be found in PARKER and 
BOWHILL (1984). 

COHERENT SCATTER DEPENDENCE ON ELECTRON DENSITY 

The radar scattering cross section is proportioned to the mean-square 
fluctuation of the retractive index n. At VHF : 

9 2 2 

n = l-Ne/emu) 
o 

where N is the electron concentration, e and m the charge and mass of the 
electron, e the permitivity of free space, and w the angular frequency of the 
impinging wave. The right-hand term is small compared to unity, fluctua- 
tions in t ) and N are proportional. Further, if the electron concentration 
increases in a scattering volume by a constant factor, the mean-square fluctua- 



tion in electron density AN^ due to turbulence will be 
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where N' is the vertical gradient of the electron density N. Thus, the scat- 
tered power P is proportional to IT within time-scales where the turbulence 
is characterized by constant mean-squared statistics. (This assumes no time- 
lag due to turbulent advection; we shall demonstrate below that this holds at 
least some of the time. 

Finally, when we divide the flare-time scattered power P at a given 
altitude by the prer- flare power P^ we obtain 

P/P = (N + AN) 2 /N 2 
o o o 

which implies 

AN/N - vfJP - 1 
o o 

D-REGION FLARE EFFECTS OBSERVED AT URBANA 

Table 1 displays features of 13 solar flare events which produced enhanced 
scatter, the enhancement shown for the altitude range 60-75 km. The sizes of 
the peak X-ray fluxes from two detectors on the GOES II satellite are shown for 
comparison. No correlation of X-ray event size and coherent-scatter power 
enhancement is evident. This is probably due to the wide variability of tur- 
bulent strength and the nonflare ion-production rate from one flare event to 
another. This suggests the need for a more sophisticated analysis. 

MODELS RELATING SOLAR X-RAYS AND ELECTRON DENSITY 

Given the GOES II X-ray measurements, we may calculate the X-ray ion-pair 
production rate at a given altitude as follows. First, we must estimate 
the X-ray spectrum from the two data points provided by the GOES detectors at 
each time. This may be done by assuming a power law spectral form: 

$ = ax b 

DO 

(see for example ROWE et al., 1970) and solving for A and B based on the wave- 
length response characteristics of the X-ray detectors (published in DONNELLY. 
1977). 


This form of the X-ray intensity ^(X) is used to calculate the desired 
ion-pair production function q according to Chapman theory. Constants 
necessary for this calculation are the average air absorption cross sections 
a (X) and the ionization efficiencies for X-rays n (X)from BANKS and KOCKERTS 
(197 3), and the scale height H and average air density M from the US STANDARD 
ATMOSPHERE (1976). 

The relationship between q (= q Q + cij and N is particularly unclear 
in the lower D region. If we assume overall charge neutrality, an unchanged 
ratio of negative ions and positive ions, and unchanged proportions of the re- 
combining species, we may derive the following (MITRA, 1974): 

XT 2 

q = aN 

where the constant of proportionality a is called the effective recombination 
coefficient (call this Case A). MITRA (1974) also proposes (Case B) 


q - 6N 
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Table 1. X-ray flare events producing measurable coherent-scatter radar 
power enhancements at Urbana between April 1978 and December 1983. 


yr 

mo 

da 

cst 

SID 

Ha 

pk flux x 10 - 5 w / m 2 
0.1~0.8nm 0.05-0.4nm 

61.5-75 km 
pk CS power 
dB above 

78 

04 

11 

0820 

3 + 

NA 

20 

3.5 

7.7 

80 

05 

21 

1510 

3 

3B 

14 

3.2 

29.8 

80 

05 

28 

135 4 

2+ 

2B 

12 

2.0 

3.5 

80 

11 

13 

1324 

2+ 

IN 

4.2 

0.7 

11.2 

80 

11 

14 

1239 

2+ 

IB 

2.4 

0.4 

13.8 

81 

01 

27 

0947 

2 

IB 

4.6 

1.3 

6.1 

81 

05 

05 

0809 

3 

3B 

12 

3.3 

2.8 

81 

08 

03 

1425 

2+ 

IN 

7.1 

1.8 

10.1 

81 

10 

14 

1111 

2+ 

IB 

30 

11 

7.4 

82 

03 

31 

1626 

2+ 

NA 

7.5 

1.7 

4.9 

82 

06 

02 

0953 

3 + 

NA 

10 

2.9 

17.4 

83 

08 

13 

1215 

1 

2B 

5.2 

0.9 

3.7 

83 

08 

21 

1159 

1 

2B 

2.2 

0.3 

6.8 


as a possible relationship, given other conditions. Both relationships require 
chemical equilibrium, which may or may not hold during a flare. If not, the 
relationship is far more complicated, but a model by THOMAS et al . (1973) pre- 
dicts delays in the lower D region between peak q and peak N of up to 14 min- 
utes, increasing with decreasing altitude (Case D) . 

DISTINGUISHING CHEMICAL SCHEME TYPES AND ESTIMATING NONFLARE ION-PAIR PRO- 
DUCTION 


To distinguish between these three possibilities (and possibly others) at 

a given altitude we may assume each possiblity in turn, compare q x at each 

time with AN/NO from the coherent-scatter data, and find the best fit. To 

simplify this comparison, note that Case A implies P«q, so that we may 

estimate q /q as 
n x H o 


q /q — P/P 
x o o 


1 


For Case B, q and N are proportional, so that 


q x /q 


o 


/p7p - l 

o 


If we make two plots, one for each of these estimates of q/q Q against the 
same values at q from the satellite data, the result would be a line of 
unit slope only for that case which is correct. If the points of the initial 
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part of the flare (ascending) do not lie along the same path as those of the 
decay phase of the flare, we may have Case D, or possibly a time lag due to 
turbulent advection. 

The result of this technique is shown in Figure 1, for the flare of Nov- 
ember 14, 1980, and the scattered power from 70.5 km. The plot of upper left 
shows a close fit to a line, but not of unit slope. The plot of upper right 
shows a fairly close fit to a line of slope 1. We may conclude q = BN, and 
more: note that the intercept of this line with the log (q / q ) = 0 line 
implieS gthjjt when q^ = q<j , log (q^) = 0.3. Thus we may estimate q = 

0.5 cm s . The lower plot shows the log scattered pcwer and ^he estimate 
of q^ based on the satellite data. The coincidence of the peaks of these 
curves validates the chemical equilibrium hypothesis, and demonstrates that 
there is no lag due to turbulent advection. 

COMBINED RESULTS FOR SEVERAL FLARES 

Figure 2 shows the type of power enhancement found for each altitude for 
each of seven flare events, arranged in order of decreasing solar zenith angle. 
The symbol E represents extremely large enhancements (some as much as 30 dB) 
which cannot be accounted for by this model. Also, note that the altitudes 
showing no response seem to proliferate downward with increasing solar zenith 
angle. This is not what Chapman theory predicts for ionization, and so pos- 
sibly indicates a variability in existence of turbulent layers. Seme events, 
particularly at high zenith angles, show delays (Case D), with delay increasing 
with decreasing altitude. Note, however, that a delay may be either a chemical 
or turbulent mixing effect. 



Figure 1. Flare-time scattered power enhancement and estimated 
electron production rate q due to flare X-rays at 70.5 km 
for the November 14, 1980 event. 




Figure 2. Types of response at each mesospheric altitude for each 
flare: A) q = aN 2 . B) q = $N. AB?) Probably A or B, but plots 
are ambiguous. D#) Delayed effect, with peak power following q 
by # minutes. E) Extreme response, power increase greater than 
cases A or B. CN) Coincident flare- time power response, but 
unable to fit to A or B model due to comparable nonflare 
fluctuations in power at that altitude. CR) Coincident response 
well above nonflare base power, but unable to fit to A or B 
models due to large fluctuations on the order of the one minute 
sample time. — ) no detectable response. 


Figure 3 shows estimates of q^ made by this technique. When the 
interpretation was unclear between Case A and Case B, but equilibrium seems to 
hold, a dashed horizontal line indicates the range of possible q estimates. 
Note that the estimates for q fall into two clusters, corresponding to small 
and large zenith angles. Also, the values of q for large zenith angle are 
smaller, as expected. Also, q Q generally increases with altitude for a given 
flare effect, also as expected. We may conclude that the estimate of q has 
some degree of reliability, although its absolute accuracy is not estimated 
here. 

DISCUSSION AND CONCLUSION 

By using models for the uncertain links between the causal X-ray emissions 
and the coherent- scatter power enhancements, we have gained insights into fea- 
tures which are not otherwise clear. Often one model fits the data better than 
others, and so we are able to choose between some features of ion-chemical 
schemes. When no existing model fits the data, as with the extreme power 
enhancements shown in Figure 3, it is unclear which part of the analysis should 
be modified: an extreme enhancement may be due to an unknown ion-chemical 

scheme, or may be due instead to sane unknown condition which allows P to 
increase faster than N 2 . This suggests further research. 




ALTITUDE (Km) 
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Figure 3, Estimates of q made by finding intercept of log (q ) vs 
log (^/q^ Plot® witfi line log C^) ^ 0 for each si titude and 
flare event which shows type A (q =aNj or type B (q =BN) 
response. Horizontal dashed lines represent the span between the 
type A estimate at q Q and type B estimate at q at altitudes 
where the choice between type A and type B is ambiguous. The 
estimate of q^ by RATNASIRI and SECHRIST (1975) is reproduced for 
comparison. 


Only altitudes with steady turbulent layers produce useful data, so 
altitude gaps are inevitable. Hcwever, the coherent- scat ter radar technique 
compares favorably with earlier methods of observing changes in N during solar 
flare, such as partial reflection and wave interaction. The temporal and 
spatial resolution are excellent for the coherent- scat ter technique, and the 
relative changes in N may be measured accurately. 
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